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Abstract. A Painleve-Gullstrand synchronization is a slicing of the space-time by 
a family of flat spacelike 3-surfaces. For spherically symmetric space-times, we show 
that a Painleve-Gullstrand synchronization only exists in the region where (dr) 2 < 1, 
r being the curvature radius of the isomctry group orbits (2-spheres). This condition 
says that the Misner-Sharp gravitational energy of these 2-spheres is not negative 
and has an intrinsic meaning in terms of the norm of the mean extrinsic curvature 
vector. It also provides an algebraic inequality involving the Weyl curvature scalar and 
the Ricci eigenvalues. We prove that the energy and momentum densities associated 
with the Weinberg complex of a Painleve-Gullstrand slice vanish in these curvature 
coordinates, and we give a new interpretation of these slices by using semi-metric 
Newtonian connections. It is also outlined that, by solving the vacuum Einstein's 
equations in a coordinate system adapted to a Painleve-Gullstrand synchronization, 
the Schwarzschild solution is directly obtained in a whole coordinate domain that 
includes the horizon and both its interior and exterior regions. 



PACS numbers: 04.20. Cv, 04.20.-q 



1. Introduction 

It is known that Painleve pQ, Gullstrand [2] and (some years later) Lemaitre [3] used a 
non-orthogonal curvature coordinate system which allows to extend the Schwarzschild 
solution inside its horizon, see Eq. ( )58|) below. In this coordinate system, from now 
on called a Painleve-Gullstrand (PG) coordinate system, the metric is not diagonal, 
but asymptotically flat, and regular across the horizon, and then, everywhere non- 
singular up to r = 0. Furthermore, one has a very simple spatial 3-geometry: the 
space-time appears foliated by a synchronization of flat instants (hereafter called PG 

| We use here the abbreviations 'PG coordinates' and 'PG synchronization' by sake of simplicity, 
without any intention of misplacing the relevant contribution by Lemaitre [3], who clarified the 



synchronizations) ||| 
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Nowadays, there is an increasing interest in the study of this type of 
synchronizations. For instance: (i) in connection with astrophysical applications, by 
taking into account that the dynamics of the gravitational collapse should be pursued 
beyond the Schwarzschild radius in a PG coordinate system [HI El [9] , (ii) in spherically 
symmetric space-times (SSSTs), as a convenient initial condition preserved under time 
evolution [TQl EI] , (hi) in relativistic hydrodynamics, when an effective Lorentzian metric 
is introduced [12], (iv) in non-relativistic situations admiting a Lorentzian description, 
namely, in 'analog gravity models' (see, for example, [131 El E3 EE EE3 IE]), and also 
(v) in modeling the black hole geometry and its associated physics, or to describe some 
quantum effects by starting from a Hamiltonian formulation [191 EH EH [22]. More 
physical issues about the use of PG coordinates and their interpretation can be found 
in [H [5l [23l [24]. For a description of the causal character of the geometric elements 
(coordinate lines, coordinate 2-surfaces and coordinate 3-surfaces) associated with PG 
coordinates see [251 12S] . 

The existence of PG synchronizations in SSSTs has been studied in the static case 
considering that the induced metric has vanishing Ricci tensor [27], and some specific 
constructions are presented in [28] . In more general cases, this existence is usually taken 
for granted but, recently, a limitation to this ansatz has been pointed out [29]. However, 
as far as the authors are aware, a definitive interpretation of this limitation as well as 
the analysis of the domains where a PG synchronization exists have not been done up to 
now. Then, some related questions arise: does every SSST admits a region of physical 
interest where a synchronization by flat instants exists? and, what are the advantages of 
adapting coordinates to a flat spatial 3-geometry? The main contribution of this paper 
is to state the above limitation clearly, in a form that it is coordinate independent, and 
to provide its physical interpretation. 

Generalized (but, in general, non-flat) PG synchronizations have been constructed 
in Schwarzschild [30] , Reissner-Nordstrom [29] , and Kerr [3T| 132] geometries, and also in 
non-vacuum SSSTs, where new insights in the study of gravitational collapse scenarios 
are achieved (by evolving an initial 3-geometry [331 El])- However, here we restrict 
ourselves to flat synchronizations in order to discuss their existence in SSSTs. 

The paper is structured as follows. Section [2] is devoted to introduce some general 
formulae for the induced geometry on space-like hypersurfaces and surfaces in SSSTs. In 
Sec. [3] the condition for the existence of a PG synchronization in SSSTs is analyzed and 
physically interpreted. In Sec. H]we write the components of the Weinberg pseudotensor 
[35] with respect to a PG synchronization, and we prove that the energy and momentum 
densities of each PG slice vanish. In Sec. [5] we consider the semi-metric connection (see 
|36j ) associated with a spherically symmetric metric expressed in PG coordinates and 
provide new insights on the Newtonian interpretation of the properties exhibited by the 
Schwarzschild field in these coordinates. Sec. [6] deals with the 3 + 1 decomposition 

coordinate character of the 'Schwarzschild singularity' obtaining an extended metric form for the 
Schwarzschild solution. In fact, the main motivation in papers [TJ [2j [3] is rather different. For historical 
remarks about this subject, and some physical interpretations see [4][5]. 
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of the Einstein equations with respect to a PG synchronization. By integration of the 
vacuum equations one recovers the extended Schwarzschild metric in PG coordinates, 
including the region inside the horizon. Finally, in Sec. [7J we discuss the role that 
our results can play for a better understanding of the geometry and physics in SSSTs. 
Some preliminary results of this work were presented at the Spanish Relativity meeting 
ERE-2009 [37]. 

Let us precise the used notation. The curvature tensor R\j of a symmetric 
connection V is defined according with the identity VjVj£ fc — VjVj£ fc = R k nj£} for 
the vector field £, and Rij = R k ikj is the Ricci tensor. We take natural units in 
which c = G = 1 and the Einstein constant is k = 8tt. We say that {t, r, 6, <p} is a 
curvature coordinate system if for constant t and r the line element is all 2 = r 2 dQ 2 
(with dfl 2 = d6 2 + sin 2 6 dip 2 the metric of the unit sphere). In these coordinates, the 
metric form is, in general, non-diagonal. When, in addition, the 3-surfaces defined 
by t = constant are flat, the curvature coordinate system is called a PG coordinate 
system. From now on, it will be understood that the title of the sections and the results 
of this work always concern space-times with spherical symmetry. The agreement for 
space-time signature is ( — h ++). 



2. Some geometrical relations 

In this section we present the geometric background needed in the following sections: 
expressions for the Ricci tensor and the extrinsic curvature of a spherically symmetric 
synchronization, as well as, the mean curvature vector and the Gauss identity for the 
2-spheres of a SSST. Of course, this material is not new and it may be bypassed or used 
as a glossary of formulae, which are conveniently referred throughout the text of the 
remaining sections. For an account on 2+2 warped space-times properties allowing to 
intrinsically characterize SSST see [38| [39]. 

Let (V4, g) be a SSST, and let us consider a canonical coordinate system {T, i£, 9, ip} 
of V4, which is adapted to the symmetries of the metric g. Then we may express the 
metric in the following general form jlQl SU H2j 

g = AdT ® dT + B dR <g> dR + C(dT ®dR + dR®dT) + D a, (1) 

with A(T, R), B(T, R), C(T, R) verifying the condition 5 = AB - C 2 < 0, D(T, R) ^ 0, 
and a = d6 <8> d6 + sin 2 9d<p> ® dip being the metric on the unit 2-sphere. 

2.1. Ricci tensor 

The induced metric 7 on the 3-surfaces T = constant, is written as 

7 = BdR ®dR + Da, (2) 

with B 7^ 0. The Ricci tensor of 7, TZic^), is given by 

( B B \ ( D F\ 

7M 7 ) = [-K-—F)dR®dR+(—K+-)(T (3) 
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where 




' - 4BD < 4) 
and TZ = 72.(7), the scalar curvature of 7, results 

dRBdRD (dnDf_^D 
ABD B r 

if <9*£>^0, (5) 
if <9 fiJ D = 0. 

Notice that the 3-surfaces T = constant are conformally nat(§] but, in general, they are 
not flat. From Eqs. (|3J), (jlj) and (J5J, we see that 7 is a flat metric if, and only if, F — 0, 
that is 

ftic( 7 ) = = (cVD) 2 . (6) 

2.2. Extrinsic curvature 

The extrinsic curvature of the slicing T = constant is defined as K = —■^£ n , -f where C n 
is the Lie derivative along the unit normal vector n, 

1 ( d C d \ _ 
n= a{dT-BdR) (?) 

with a 2 = — — . From (j2J), we obtain 



where 



K = ^BdR®dR + <&Da (8) 



* = = JL f 20*0 - - d T B ) (9) 



$ = K e e = K« = — ( —OrD - d T D ) (10) 



are the eigenvalues of i^. Developing the Lie derivative of K with respect to the shift 

Q_d_ 

B dR 



vector (5 = we arrive to the expression 



C p K=^d R 



c ,2 



mB 



c 

dR®dR + —d e (<$>D)o-, (11) 
B 



B 

which will be needed in Sec. [6] to split the Einstein evolution equations with respect to 
a PG synchronization. 

§ Notice that the Cotton tensor of 7, Cijk{^) = DiQjk — DjQn~, {Di is the covariant derivative with the 
Levi-Civita connection of 7, and Qij = TZij — ^"fij) identically vanishes, djk(j) — 0. This is the result 
that one could expect ought to the algebraic properties of the Cotton tensor and the assumed spherical 
symmetry. This means that 7 is a conformally flat metric. Then, 7 may always be written in isotropic 
form by making a coordinate change on each 3-surface T = constant. An interesting summary on the 
Cotton tensor properties is given in [43]. 
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2.3. Mean curvature vector 

The mean curvature vector if of a 2-sphere S defined by constant T and R is given by 

H = ~D6 Y B9tD ~ CdRD) fff + {A9rD ~ CdTD) dRl ' (12) 
This expression directly follows by taking the trace (with respect to the induced metric 
Da) of the extrinsic curvature tensor K, of each S, which is defined by 

/C(e a , e b ) = - (V eo e 6 ) ± = - (v T ab ± + I* A) (13) 

where the minus sign is taken as a matter of convention. V is the Levi-Civita connection 
of g, {ei = is a coordinate basis of S, i = a,b = 9, ip, and the symbol - 1 stands for 
the projection on S ± : the time-like 2-surface orthogonal to S. For a detailed study of 
/C with applications in relativity see, for example, [3HH21HS]- Then, it results 

K = l -Do ® H. (14) 
The one form T metrically equivalent to H, T a = g a pH@, is written as 

T = -—(d T DdT + d R DdR) = -dhxD, (15) 
which can be also obtained from the general expressions presented in |4T] . 



2.4- Gauss identity 

Given a space-like 2-surface S of a space-time (V4, g), the Gauss identity provides a scalar 
relation involving the background geometry and the intrinsic and extrinsic properties of 
S, and it may be expressed adjj] 

R(h) = -R(g) + g(H, H) + 2tr(i^ x K k ) + 2Ric(l, k) - 2W(l, k, I, k) (16) 
3 

where R(h) is the scalar curvature of the induced metric h on E, R(g), Ric and W are, 
respectively, the scalar curvature, the Ricci and the Weyl tensor of g, and I and k are 
two independent future pointing null vectors normal to S satisfying g(l, k) = —1. Now, 
H is the mean curvature of E, that is H = trJC, with /C the extrinsic curvature tensor 
of E (as defined by the first equality of ( 1TB1 . taking e a and e& tangent to E); Ki and 
Kk are, respectively, the second fundamental forms of E relative to I and k, that is, 
Ki(e a ,eb) = g(l,K,(e a ,eb)), and similarly for Kk. The trace is taken with respect to h, 
that is, tr(K t x K k ) = h^h^K^K^. 

|| The Gauss identity is usually given in terms of the Riemann curvature tensor R a /3-yS (see, for instance 
[45l|46]) from which the expression ((T6|) follows by taking into account the algebraic decomposition: 

R 

RafljS = W a fSjS + 9a[jRs]P — 9/3[-yRs]a ~ ~^ 9a[~j9S}[3 

where the bracket stands for antisymmetrization of index pairs, T[ a p] = \ {T a p — Tp a ). 
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For £ — S, a 2-sphere (orbit) in a SSST, the Gauss identity may be written as 

p = 2\+ l -{Ap-p l -p 2 )+ l -H 2 . (17) 

with p = 2/D the scalar curvature of Da, A the simple eigenvalue^! of W (provided that 
W 7^ 0), and {^,^1,^2} stand for the Ricci eigenvalues. More precisely, the tangent 
space to S is the eigenplane associated with the eigenvalue p. The eigenvalues p\ 
and P2 may be real (simple or multiple) or complex, and have associated the invariant 
2-plane tangent to S 1 - . Here, we have considered the algebraic classification of Churchill- 
Plebahski for a symmetric 2-tensor in Lorentzian geometry, and its peculiarities in 
spherical symmetry (see jlHj, and also jl9] for an intrinsic approach). We have taken 
into account that R(g) = 2p + pi + p 2 , and 2Ric(l,k) = —(pi + p 2 )- In addition 
H 2 = g(H,H) = —4tr(Ki x K^) because, according to ( fl4l) . the second fundamental 
forms Ki and Kk are both proportional to h = Da. 



3. Painleve-Gullstrand slicings 

In this section we find the condition ensuring the existence of a PG-synchronization in 
spherical symmetry, and we discuss its invariant meaning in terms of the eigenvalues 
of the Weyl and Ricci tensors. Also, using a radial curvature coordinate r, we provide 
a physical interpretation of this condition in terms of the Misner-Sharp gravitational 
energy of a 2-sphere of radius r. As commented below, this result shows an interesting 
interconnection between the study and classification of trapped surfaces (see [50]) and 
the existence of flat slicings. 



3.1. Existence condition 

Let us start from the general metric form (CQ). Exploring the gauge freedom to make 
coordinate transformations of the form T = T(t,r), R = R(t,r), we look for a function 
t(T,R) whose level hypersurfaces t = constant are Euclidean, i. e., the induced metric 
is positive and flat. Under such a transformation, the metric Q is expressed as 

ds 2 = £ 2 dt 2 + X 2 dr 2 + 2£- X dtdr + V(t, r)dtt 2 , (18) 

with V(t,r) = D(T(t,r), R(t,r)), and the vector fields £ and \ defined by 

Over-dot and prime stand for partial derivative with respect t and r, respectively, and 
J = TR' — T'R ^ assures coordinate regularity. 

The scalar products £ 2 = g((,,£,), x 2 = fi'(X) X) an d £ " X = x) can b e written as 

5(dD) 2 f 2 + 2VZT + BV 2 = (d R D) 2 £ 2 , (20) 

If One has 

\w^{l a k p ~ fk a ) = Xil^K - l u kp), 
and then, W(l, k, I, k) = W^^l^l 01 ^ = -A. 
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S(dD) 2 T' 2 + 2D' ZT' + BV' 2 = {d R D) 2 x 2 , (21) 

(d R D) 2 £ ■ X = 5(dD) 2 TT' + Z (V'f + VT') + BVV 1 , (22) 

where we have used the relations D = Td^D + Rd R D and V = T'drD + R'd R D to 
substitute R and R' in terms of T and T' . We have denoted Z = Cd R D — B8tD. Then, 
due to the Lorentzian character of the metric (5 < 0), Eqs. ( l20l) and (l2~TI) lead to real 
values for T and T' if, and only if, the inequalities 

i 2 (dV) 2 < V 2 , x 2 (dV) 2 < V' 2 (23) 

are satisfied. Now, looking for a flat synchronization, we have that the induced metric 
on the 3-surfaces t = constant is flat if, and only if, 

4V X 2 =T>' 2 (24) 

according to Eq. ([6]). Consequently, in the case of a flat synchronization, the second 
inequality in ( 1231 is equivalent to 

(dVv) 2 < 1. (25) 

So, under the assumed spherical symmetry, Eq. (125]) provides the necessary and 
sufficient condition to be fulfilled for the existence of a flat slicing. The first inequality 
in (|23|) guarantees that the slices are space-like, that is, that the slicing is a PG 
synchronizat ion . 

3.2. Geometric interpretation 

In terms of the scalar curvature p = 2/T> of the metric T>dQ 2 , the above condition (T25l) 
may be expressed as follows 

(dp) 2 < 2p 3 , (26) 



which involves the sole invariant p. On the other hand, according to Eq. (T15]) . 
H 2 = T 2 = (dlnV) 2 , and then Eq. ( 125]) gives an upper bound for the norm of the 
mean extrinsic curvature H of the group orbits, 

R2 <^ = 2P- (27) 
Moreover, from the Gauss relation (fITj) . we arrive to the following result. 

Proposition 1 In a spherically symmetric space-time the following conditions are 
equivalent. 

(i) There exists a Painleve-Gullstrand synchronization. 

(ii) (dp) 2 < 2p 3 , where p is the scalar curvature of the 2-spheres. 

(Hi) H 2 < 2p, where H is the mean curvature vector of the 2-spheres. 

(iv) pi + /i 2 — 4/i < 6A, where p.\,p%, and p (double) are the Ricci eigenvalues, and 
A is the simple eigenvalue of the Weyl tensor, or A = when the space-time is 
conformally fiat. 
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Notice that this is a geometric result, which will be physically interpreted in the next 
subsection. Taking into account the Einstein equations, in the above item (iv) the 
Ricci eigenvalues, {^x, ^2, H>} may by substituted by the corresponding energy tensor 
eigenvalues, {ei,e 2 ,e}, giving 

ei + e 2 - e < 3A. (28) 
3.3. Physical interpretation 

By definition, see [HI r is a coordinate of curvature for the spherically symmetric 
metric form (flHjl if V(t,r) = r 2 , so that V = 2r and T> — 0. Then, (1251) says that 
((ir) 2 < 1, and taking into account that the Misner-Sharp gravitational energy E of a 
2-sphere of radius r is expressed as (see [51], [50] ) 

£=^(l-(dr) 2 ), (29) 

we arrive to the following result. 

Proposition 2 Any spherically symmetric space-time admits a Painleve-Gullstrand 
synchronization in the region where the Misner-Sharp gravitational energy is non- 
negative, E > 0. 

The Misner-Sharp energy has been painstakingly analyzed in [50J, providing useful 
criteria to study trapped surfaces. The main novelty here has been to relate this concept 
and the existence of PG synchronizations^ 

Moreover, the flatness condition (pHl) implies that x 2 = 1, and the metric (TT8l) is 
written as 

ds 2 = £ 2 dt 2 + 2£-xdtdr + dr 2 + r 2 dtl 2 . (30) 
Then, accordingly to f[2"0"l) and fT22l . the following relations must occur 
A(t,r) =e = J 2 5(dr) 2 

(31) 

B(t,r)^^x = eJ^6[(dr) 2 -l} 

where e = ±1. So, the real function B exists in the region where (dr) 2 < 1, and we have 
the following result. 

Proposition 3 Let r be the radius of curvature of the orbits (2-spheres) of the isometry 
group of a spherically symmetric space-time with metric g. In the region defined by the 
condition 

(dr) 2 = g^d^r < 1, (32) 

the Misner-Sharp energy is not negative and a curvature coordinate system {t, r, 8, (p} 
exists in which the metric line element may be written as 

ds 2 = A(t, r) dt 2 + 2B(t, r) dt dr + dr 2 + r 2 dQ 2 , (33) 

where dfl 2 = dQ 2 + sin 2 9 dip 2 . 

+ Marc Mars inspired us in the obtaining of this relation. 
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The Lorentzian character of the metric impose that the functions A and B must 
satisfy the sole restriction A < B 2 , which is implied by fl3Tj) . 

The Misner-Sharp energy is a geometric invariant that may be physically interpreted 
as an effective gravitational energy whose origin is the interaction between the energetic 
content and its associated field (see [50]). Given that the intrinsic and extrinsic scalar 
curvatures of the 2-spheres are p = 2/r 2 and H 2 = 4(<ilnr) 2 , and according with (1291) 
one has the invariant expression 

£=-L(l-— H 2 ). (34) 

Finally, notice that Eq. ( 1271) does not constraint the causal character of the mean 
curvature vector H, which might be time-like, light-like or space-like. This is a 
remarkable property, because a 2-sphere is said to be trapped, marginal or untrapped if 
H is, respectively, time-like, light-like or space-like (see e. g. [4"5 | |4"6 | |50 | 152]). 



4. Energy and momenta densities of a Painleve-Gullstrand slice 

In this section we establish the following result. 

Proposition 4 In any spherically symmetric space-time, the Weinberg energy and 
momenta densities vanish for every Painleve-Gullstrand synchronization. 

Of course, to find a coordinate system in which the Weinberg densities vanish is 
not a surprising property, due to the non-tensorial character of them. However, the 
novelty here is to show that, for every SSST, such a vanishing property occurs in PG 
coordinates. 

In order to proof the above result, let us consider the metric ( I3B1 written in a 
quasi-Minkowskian form, that is g = rj + h with rj the Minkowski metric, h 00 = 1 + A, 
hoi = Bxi/r, and hij = 0. 

We start from the expression of the Weinberg pseudo-tensor [35] , 

2Q i0X = _ %oa _ dhP* + dh^ v ox + dh^_dh^^ 

dxo dxi dx^ dx^ dxi 8xq ' 

where Latin and Greek indexes go from 1 to 3 and from to 3, respectively, and all 
indexes are raised and lowered with the flat metric t]. In this case, it results Q l0 ° = 
(according with [53]) and 

2Q"» =(^ + ff)s ij +(-- B') ^ (35) 
\r J \r J r r 

The derivative of this expression leads to 



dQW ( B' B „.A .xi. (B B' 
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and by contraction of the indexes, directly follows that = 0. Then, the four- 

momentum density vanishes, 

^ 0A = = (36) 

and hence, the angular momentum densities ] %x = x l r ox — x x r° l also vanish, according 
with the announced conclusion. 

For the special case of the Schwarzschild geometry, the vanishing of the energy 
density may be intuitively understood invoking the Einstein equivalence principle. 
Taking e = 1 in the extended form ( 158]) of the Schwarzschild metric, t represents the 
proper time of a radial geodesic observer which initially stays, in r = oo, at rest with 
respect to a static observer. Locally, such an observer does not feel any gravitational 
effect. 

5. Painleve-Gullstrand slicings and semi-metric connections 

In the eighties, Bel proposed an extended Newtonian theory of gravitation based on 
a semi-metric connection associated with an observer congruence and a flat spatial 
3-metric [36]. In a space-time, with metric g^, which admits a spatially flat slicing 
given by the coordinate hypersurfaces x° = constant, the connection coefficients of the 
aforementioned semi-metric connection are written as [36] , 

A fc = - Ito = U kl {d l9m - 2d o9oi ) (37) 

ft* = - 21* = 5 ki (d igoj - 2d jgoi ). (38) 

Consequently, a SSST metric admits a Newtonian interpretation when it is written 
in PG coordinates and it is considered in the above context. In fact, taking into account 
the expression fl33|) of the metric, we have g 00 = A, goi = Bxi/r, and then 



d ig0 j = —S-ij + ( B' ) — ^ = djg 0i 



Then, the connection coefficients result 



A k = -(a'-2B) — (39) 



2V J r 

ftf = 0, (40) 

which means that, in the region of a SSST where a PG synchronization exist, the 
gravitational field may be interpreted as an inertial field of radial accelerations and 
vanishing rotation. 

In particular, for the case of the Schwarzschild metric, we have A = — (l — ^r), 
B = , and then the vector component of the connection reduces to 
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where e r is the unit vector in the radial direction. The above expression (1411) gives the 
acceleration of a unit mass particle radially falling in the Newtonian field of a mass m. 
Similar Newtonian interpretations have been considered from a different point of view 
(see, for example, [U [5j [23]). 



6. Painleve-Gullstrand slicings and Einstein equations 

In General Relativity, when dealing with the evolution (or 3 + 1) formalism (see [54], [55]. 
and [56] for a recent review) one introduces a vorticity free observer n, n 2 = —1, and 
Einstein equations are decomposed in the following set of constraint equations (k is the 
Einstein constant), 

Tlia) + (trA^) 2 - trA^ 2 = 2kt (42) 

V- {K-trKj) = Kq (43) 
and this other system of evolution equations 

da = -2aK + £^7 (44) 

8 t K = -Wa-Ka[U+-(T-p)j}+a[TZic(j)+tTK K-2K 2 ]+C P K . (45) 

Here, 7 and K are, respectively, the metric and the extrinsic curvature of the space- 
like slices whose normal vector is n; V is the Levi-Civita connection of 7, and the Ricci 
tensor and scalar curvature of 7 are denoted by lZic{a) and ^(7), respectively; the trace 
operator associated with 7 is denoted by tr, so that, (V • K) a = (trVA") a = j 13 "^iKja is 
the divergence of K with respect to 7. In the usual evolution formalism notation, n is 
written asn = a_1 (^ — where a is the lapse function and (3 is the shift vector. 
The energy content T = {r, q,p, 11} has been decomposed relatively to n, that is 

T=rn®n + n®q + q®n + Yl + pj, (46) 

with r = Tin, n), q = — _l_T(n, •), p and II being the energy density, the energy flux, the 
mean pressure and the traceless anisotropic pressure as measured by n, respectively; _L is 
the projector on the 3-space orthogonal to n associated with the 3-metric 7 = g + n®n. 



6.1. Spherical symmetry 

In the case of a SSST, using the expression (jSJ) of the extrinsic curvature, the constraint 
equations (}4"2"]) and (14"3]1 are equivalent to 

$($ + 2^) = KT - y (47) 

2d R 3> + - = - Kq R (48) 
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where q R is now the radial component of the energy flux. For the evolution equation 
5}, taking into account the expression (jlip . we have IT W = Hqq sin 2 6 and 



&r(*B) = -VBd R (^\ -Ka(u RR +^(r-p)B 



5 



5* 



(49) 



+a (jK + ^ + Z>) + -d R {D$). 



(50) 



For the metric (CD), Eqs. (I7J), IjiBjl . (139")) and (150]) are the 3 + 1 splitting of the Einstein 
equations with respect to a vorticity free observer. The proper space of such an observer 
is Euclidean if, and only if, F = 0, and then 1Z = 0. The integration of these 
equations for simple energetic contents (for instance, a dust model) should provide the 
corresponding metric form in PG coordinates. In the next section, the vacuum case is 
considered: the extended form of Schwarzschild solution in PG coordinates is obtained 
from the sole consideration of the field equations. 

6.2. Schwarzschild vacuum solution 

The extended Painleve-Gullstrand-Lemaitre metric form of the Schwarzschild solution 
may be obtained assuming spherical symmetry and the existence of a flat 
synchronization, TZic( , ~f) = 0, and then, solving the vacuum Einstein equations in a 
coordinate system adapted to such a synchronization. So, let us take r = p = Hrr = 
Hgg = F = 1Z = 0. Then, for the metric expression ( 133]) . the lapse function is given by 
a 2 = B 2 — A and the constraint equations (1471) and (148]) result 

$($ + 2$) =0 (51) 
r $' + $ _ $ = o (52) 

with 

IB B' 

$ = --, * = — . 53 

a r a 

When $ = 0, taking into account also the evolution equations, we recover the 
Minkowski space-time. In the generic case, <3> = —2^> ^ 0, Eq. (153]) leads to 

B = f(t)r~ 1 / 2 (54) 

with fit) an arbitrary function. Substituting Eq. ( 154]) in the momentum constraint 
(152]) . it reduces to a' = 0. Consequently, the lapse is a function of the sole variable t, 
a(t), and we can take a = 1 by re-scaling the coordinate t. Then, we have 

$ = f(t)r~ 3/2 = -2^. (55) 
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Next, the evolution equations (I49p and (1501) are written as 



^ = ^(2$ - tf) + ^(B 2 ^)' 



$ = + — (r 2 $)'. 



r 



Given that $ 



— 2^, these last equations are equivalent to 



r 



(56) 



3# + B' = 0. 



(57) 



r 



By using the expressions ( 1541) and ( 1551) . the equation ( 1561) leads to f(t) = constant and 
the equation (1571) is identically satisfied. Finally, by taking / = e\/2m, we obtain 



which is the extended form of the Schwarzschild solution obtained by Painleve, 
Gullstrand and Lemaitre [B 12 [3]. The positive parameter m is the Schwarzschild 
energy. The sign e provides two coordinate branches for the solution, the Kruskal- 
Szekeres black and white hole regions being described by the above metric with e = 1 
and e = — 1, respectively, see [21]. Note that the r coordinate can take any positive 
value, < r < +oo. In fact, from ( 158]) . we have [air) 2 = g rr — 1 — — < 1, and the 
domain of a Painleve-Gullstrand chart extends for every value of r ^ 0. Notice that 
(T2"9~j) implies that E = m, which provides the physical interpretation of the parameter 

m as an effective energy [50]. Moreover, writing dt = dts + ey^p(l — 2 y 1 )~ 1 dr, one 
recovers the usual Schwarzschild metric form 



where t$ is the coordinate time of the static observer (— oo < t$ < oo) and the rank of the 
r coordinate is restricted to be r > 2m. According to the Jebsen-Birkhoff theorem (see 
[57] an references therein), we recover the Schwarzschild metric as the sole spherically 
symmetric solution of the vacuum Einstein equations. 

Other derivations of the Schwarzschild solution providing improvements of the 
original proof of the Jebsen-Birkhoff theorem have been achieved by solving the 
field equations in null coordinates (see j50j EH] and reference therein). From the 
conceptual point of view, any of these derivations make unnecessary to get a coordinate 
transformation allowing to extend the domain of Schwarzschild chart from the outer to 
the inner horizon regions. 




(58) 




(59) 
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7. Discussion 

In this work we have analyzed the existence of flat synchronizations in SSSTs. The 
condition fl27j) provides an upper bound for the norm of the mean extrinsic curvature 
vector of the isometry group orbits which, using the Gauss identity, may be expressed 
in terms of curvature invariants. Moreover, the associated flat slices have vanishing 
Weinberg energy and momentum densities. We have seen that any spherically symmetric 
metric admits a Newtonian interpretation in the context of the Bel extended Newtonian 
theory of gravitation. Our study offers a new perspective about the meaning of the 
Painleve-Gullstrand coordinates. This study applies for any SSST in the region where 
these coordinates exist. In this region, the gradient of the radial Painleve-Gullstrand 
coordinate r may be space-like, light-like, or time-like, according to the condition 
(dr) 2 < 1, which means that the Misner-Sharp gravitational energy of a sphere of radius 
r is non-negative. This condition may be tested for any SSST, starting from the general 
metric form (CQ). For instance, it occurs elsewhere in the Schwarzschild geometry, as it 
has been pointed out at the end of Sec. [6j Moreover, one has that (dr) 2 < 1 everywhere 
for any Robertson- Walker metric with an energetic content which satisfies the usual 
energy conditions. In fact, if we put in (GD) A = -1, B = a 2 (t)/(l + \r 2 ) 2 (with 
k — 1, 0, — 1 the universe curvature index), C = and D = r 2 B, we obtain that, in this 
case, ( J27|) is equivalent to k + a 2 > 0, which means that the proper energy density of 
the cosmological fluid is non negative. Consequently, any Robertson- Walker space-time 
that satisfies this energy condition admits a PG synchronization. This property is also 
obtained directly from the inequality (128]) . In this case, r = —e\ and p = ei = e are, 
respectively, the energy density and the pressure of the cosmological fluid, and A = 0, 
because the Robertson- Walker metric is conformally flat. We leave for a future work 
the obtaining of the Painleve-Gullstrand form of these Robertson- Walker cosmological 
models. 

Finally, we have presented an improved proof of the Jebsen-Birkhoff theorem by 
expressing and solving the vacuum Einstein Equations in PG-coordinates. So, the 
extended Painleve-Gullstrand-Lemaitre metric form of the Schwarzschild solution is 
directly obtained. 
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